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Systemzystairodleini 

Die Erfindung betrifft ein Sensorsystem und ein Verfahren zur Ermittlung von 
Systemzustanden, im folgenden Zustande genannt. 

Bei Systemen, deren Zustandsgrossen mittels zumindest eines Sensors eines 
Sensorsystems mit einer zugeordneten Recheneinheit zu erfassen sind, ist die 
moglichst genaue Bestimmung derselben aus zumindest einer MeBgroBe je Zeitpunkt 
von groBer Bedeutung fur die Funktionsfahigkeit des Systems. Dies gilt insbesondere 
fur sicherheitskritische Systeme, beispielsweise Luftdatensysteme eines Flugzeugs, 
da bei nicht ausreichender Genauigkeit der MeBgrossen weitere MaBnahmen, wie z.B. 
weitere Sensoren, erforderlich sind, um die geforderte Sicherheit des Gesamtsystems 
zu erreichen. Weitere Beispiele finden sich bei der Bestimmung von Zustanden bei 
chemischen und nuklearen Reaktionen oder in der Navigation bei der Bestimmung der 
Position aus mehreren Steckenmessungen wie GPS oder Radar. 

Bei Systemen mit Sensorsystem mit zumindest einem Sensor und diesem 
zugeordneten Recheneinheit zur Ermittlung zumindest einer ZustandsgroBe des 
Systems sind Verfahren bekannt, die iiber die Invertierung einer eindimensionalen 
Kalibrierkurve eines Sensors dem MeBwert des Sensors gegebenenfalls in 
Abhangigkeit weiterer Parameter des Systems eine bestimmte ZustandsgroBe 
zuordnet. Beispielsweise erfolgt bei einem Luftdatensystem eines Flugzeugs die 
Messung des momentanen Flugzeug-Zustands mittels Luftdaten-Sensoren, die z.B. 
statische und dynamische Driicke in der Nahe des Flugzeug-Rumpfes messen. Uber 
die eindimensionale Invertierung der zuvor bestimmten Kalibrierungen der Messung in 
Abhangigkeit eines Zustandes, also z.B. eine Druckmessung in Abhangigkeit des 
statischen Druckes, werden aus diesen MeBgrossen Zustandsgrossen der ungestorten 
Stromung, wie z.B. Stromungswinkel, Machzahl, Hohe und dynamischer Druck, 
ermittelt. 

Nachteilig an diesem Verfahren ist, daB die eindimensionale Invertierung von 
mehrdimensionalen Kalibrierungen nur eine ungenaue Ermittlung der Zustande zulaBt. 
Auch konnen Stabilitatsprobleme auftreten. 

Es wiirde deshalb ein Verfahren fur ein Luftdatensystem entwickelt, bei dem mit Hilfe 
einer Kostenfunktion x 2 (x> y, u) , die Kalibrier-Kurven oder -Flachen umfaBt, eine 
Bestimmung der momentanten Zustandsgrossen x des Systems aus zumindest einem 
Sensorsignal vorgenommen wird, um die Genauigkeit und Zuverlassigkeit der zu 
bestimmenden Zustande x zu optimieren. Dieses Verfahren ist beispielsweise in 
Schatt, 18.07.2000 
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(Friehmelt H., Jost M., Flush Air Data System- Fortschrittliches Luftdatensystem fur 
die Luft-und Raurnfahrt, DGLR Jahrestagung 1999, Berlin, Vortrag 99-180, Seite 5) 
veroffentlicht. Bei diesem Verfahren werden Kalibrierkurven Oder -flachen des 
Sensorsignals y, der allgemein ein Vektor ist, und von den gesuchten 
Zustandsgrossen x und bekannten Konfigurationen bzw. Reglereingaben u abhangt 
verwendet 

Bei dem Luftdatensystem, fur das dieses Verfahren vorgeschlagen ist, stellt x die 
Luftdaten a, p, q c und p s , y die Druckmessungen, d.h. die MeBgrossen, und u 
Flugzeugkonfigurationen, also systemabhangige Paremeter, dar. 

Zur Bestimmung der momentanten Zustandsgrossen x ist in dem vorgeschlagenen 
Luftdatensystem eine Minimierung der Kostenfunktion % 2 (x, y, u) mittels eines 
Gradientenabstiegsverfahrens vorgesehen. Die Kostenfunktion x 2 ( x o, y, u) wird 
ausgehend von einer zufallig gewahlten Anfangsinitialisierung x 0 des Zustandes durch 
eine rekursive Variation von x entlang des Gradienten von % 2 so lange erniedrigt bis 
ein lokales Minimum von x 2 bezuglich x erreicht wird. 

Nachteilig an diesem Verfahren ist, daB mehrere Iterationszyklen zum Erreichen des 
globalen Minimums notwendig sind, wodurch die Rechenzeit des Systems erhoht 
wird. Auch ist nachteilig, daB das Erreichen eines globalen Minimums von % 2 nicht 
garantiert werden kann. Das beschriebene Verfahren ist insbesondere in solchen 
Fallen nachteilig, bei denen die Kalibrierflachen sehr zerkluftet sein konnen. Dies tritt 
bei einem Luftdatensystem z.B. in der Transonik auf. 

Es ist daher die Aufgabe der Erfindung, ein Verfahren zur Implementierung in einem 
MeBsystem sowie ein MeBsystem bereitzustellen, welches die Bestimmung des 
Zustandswertes aus einem MeBwert weiter optimiert und garantiert. 

Diese Aufgabe wird mit den Merkmalen der unabhangigen Anspruche gelost. Weitere 
Ausfuhrungsformen sind in den Unteranspruchen angegeben. 

Im fblgenden wird die Erfindung anhand der beiliegenden Figuren beschrieben, die 
zeigen: 

Figur 1 eine schematische Darstellung eines Systems mit einem Sensorsystem, 
fur das das erfindungsgemaBe Verfahren vorgesehen ist, 

Figur 2 eine Darstellung einer Abfolge des erfindungsgemaBen Verfahrens. 

Das in der Figur 1 dargestellte System 1 weist ein Sensorsystem 2 mit zumindest 
einem Sensor 3 und einer Recheneinheit 4 auf. Der Sensor 3 liefert zumindest einen 
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MeBwert, aus dem die Recheneinheit 4 eine Oder mehrere Zustandsgrossen ermittelt. 
In der Recheneinheit 4 ist das erfindungsgemaBe Verfahren implementiert und kann 
zusatzlich eine Regelungsvorrichtung zur Regelung von Zustanden bzw. 
Zustandsgrossen des Systems 1 umfassen. Die Recheneinheit 4 erstellt Stellsignale 
oder Kommandos fur eine Stelleinheit 7 zur Veranderung zumindest einer 
Zustandsgrosse des Systems 1. Eine Anderung des Systemzustands, die mit dem 
Bezugszeichen 10 bezeichnet ist, wird wiederum von dem zumindest einen Sensor 3 
erfaBt. 

Das System 1 kann fur eine Vielzahl von Anwendungen vorgesehen sein, z.B. die 
Ermittlung von Zustandsgrossen eines Flugzeugs aus Druckmessungen, die 
Bestimmung von Konzentrationen oder Temperaturen bei chemischen Prozessen, die 
Ortsfindung in der Navigation durch Vergleich verschiedener GPS oder Radar 
Abstandsmessungen oder die Bestimmung von Deformationen und Kraften in der 
Strukturmechanik. 

In der Recheneinheit 4 sind zur Durchfuhrung des erfindungsgemaBen Verfahrens 
Kalibrier-Flachen oder -Kurven zur Verwendung in einer Kostenfunktion abgespeichert, 
die mit Hilfe von Versuchen ermittelt worden sind und die einen funktionalen 
Zusammenhang darstellt, mit dem aus zumindest einer ZustandsgroBe MeB- oder 
Sensor-Werte bestimmt werden. 

Die Kostenfunktion x 2 (x> y> u) wird in einem Verfahrensschritt 10 vorzugsweise in der 
Recheneinheit 4 in einer Weise gebildet, so daB sie die Abweichung eines 
tatsachlichen MeBwerts von den Kalibrierkurven oder -flachen, im folgenden kurz 
Kalibrierungen genannt, in Abhangigkeit des unbekannten Zustands angibt. Dabei 
kann die Kostenfunktion eine Gewichtung einer Abweichung des MeBwerts von der 
Kalibrierung durch die MeBwert- bzw. Kalibriergenauigkeit enthalten. Die verwendeten 
GroBen x, y und u sind vorzugsweise vektorielle GroBen. 

Fur ein Luftdatensystems eines Flugzeugs kann die Kostenfunktion beispielsweise in 
folgender Form vorliegen: 

X 2 = [y-(c P q c + Ps)] T [Cov(Cp)q c 2 +a p 2 r 1 [y-(c p q c +p s )] 

Dabei sind y die Druckmessungen, c p (x, u) zuvor bestimmte Kalibrier-Kurven oder - 
Flachen in Form eines funktionalen Zusammenhangs und Cov(Cp)(x, u) zuvor 
bestimmte Genauigkeiten der Kalibrierflachen. Der Zustand x wird aus den Luftdaten 
a, p, q c und p s gebildet und a p beschreibt die Druckwandlergenauigkeit der Druck- 
Offnungen. 
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Fur das erfindungsgemaBe Verfahren fur jegliche Anwendung ist ein entsprechender 
funktionaler Zusammenhang zwischen MeBgroBen und Systemzustanden in einer 
Speichereinheit abgelegt und wird zusammen mit den Systemparametern u in einem 
Schritt 13 zur Bildung der Kostenfunktion der Recheneinheit 14 zugefuhrt. Weiterhin 
wird der Recheneinheit die jeweilige Messung 17 zugefuhrt. Diese Kostenfunktion ist 
fur jede Messung y, also z.B. in jedem Zeitschritt, neu zu erstellen, da sie abhangt von 
den jeweiligen MeBwerten y und gegebenenfalls dem jeweils relevanten System- 
Parameter u. 

Fur jede Kostenfunktion einer Messung wird das Minimum der Kostenfunktion wie 
folgt ermittelt: 

Vorzugsweise ist vor Ablauf des erfindungsgemaBenVerfahrens bzw. im System 1 eine 
Approximationsfunktion der Kostenfunktion festgelegt, die so beschaffen ist, daB 
Minima derselben ermittelt werden konnen. Alternativ konnen wahrend des 
erfindungsgemaBen Verfahrens auch mehrere Approximationsfunktionen in 
verschiedenen Zustandsbereichen zur Auswahl stehen, die in einem Schritt 19 
ausgewahlt werden. In einem Schritt 20 wird die Approximationsfunktion auf der Basis 
der festgelegten Funktionsart fur vorzugsweise vor Ablauf des Verfahrens festgelegte 
Approximationsbereiche aufgestellt. Die Festlegung der Approximatonsbereiche hiangt 
von der Kostenfunktion und der ausgewahlten Approximationsfunktion sowie einem 
maximal tolerierbaren Fehler iiber den gesamten relevanten Zustandsbereich ab. Die 
Summe der Approximations-Bereiche decken dabei den gesamten relevanten 
Zustandsbereich ab. Dabei kann auch nur ein Approximationsbereich fur den 
gesamten relevanten Zustandsbereich vorgesehen sein. 

Wenn alternativ die Festlegung der Approximations-Bereiche im Verlauf des 
erfindungsgemaBen Verfahrens bzw. im System 1 erfolgt, werden die Approximations- 
Bereiche innerhalb des Zustandsbereichs x in einem Schritt 30 angepaBt. Dies 
geschieht mit rekursiven Verfahren, wobei beispielsweise eine weitere zur Verfugung 
stehende Approximationsfunktion verwendet wird oder der Approximations-Bereich 
reduziert wird. Dabei konnen die Ansatze der bekannten Verfahren der neuronalen 
Netze zugrunde gelegt werden. 

Die Anzahl der bei der Minimum-Bestimmung verwendeten Approximations-Bereiche 
kann auch bei der Verwendung von speziellen, z.B. quadratischen 
Approximationsfunktionen reduziert werden, wenn ausgeschlossen werden kann, daB 
in einigen Approximations-Bereichen kein Minimum vorliegen kann. 

In alien Approximations-Bereichen, in denen sich Minima befinden konnen, werden 
dann in einem Schritt 40 alle Minima in diesen Approximations-Bereichen bestimmt. 
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Sofern dabei im relevanten Zustandsbereich mehrere lokale Minima ermittelt worden 
sind, erfolgt in einem Schritt 50 eine Bestimmung des globalen Minimums aus dem 
Vergleich der lokalen Minima in einer Zustands-Bestimmung 25. Dabei wird durch 
einen Vergleich der Kostenfunktionen % 2 W aller lokalen Minima x uber alle 
Approximations-Bereiche das globale Minimum von x 2 (*) identifiziert. Die 
nachfolgende Zustands-Bestimmung 60 ermittelt dann die dem jeweiligen MeBwert 
zugehorige Zustandsgrosse x. Im Beispiel der Lufdatenmessung ist das Ergebnis 
dieser Minimierung der gesuchte Luftdatenvektor x mit den Luftdaten a, p, q c und p s . 

Vorzugsweise liegt bei dem erfindungsgemaBen Verfahren die Approximationsfunktion 
sowie die Einteilung der Kostenfunktion in Approximations-Bereiche fest. 
Entsprechende funktionale Vorkehrungen konnen in der Recheneinheit 4 zur 
Verfiigung stehen. 

Wenn bekannt ist, daB eine spezielle verwendete Approximationsfunktion in einem 
Approximations-Bereich nur ein einziges Minimum haben kann, ist zur Ermittlung des 
lokalen Minimums in einem Approximations-Bereich ein Gradientenverfahren mit Hilfe 
eines beliebigen oder aufgrund einer Abschatzung festgelegteri Startvektors im 
jeweiligen Approximations-Bereich vorteilhaft. Alternativ kann auch ein 
Halbierungsverfahren oder ein vergleichbares Verfahren angewendet werden. 

Bei entsprechenden Approximations-Funktionen, bei denen auch mehrere Minima in 
dem relevanten Approximations-Bereich des betreffenden Gitterpunktes vorliegen 
konnen, konnen erfindungsgemaB analytische Verfahren zur Minimum-Bestimmung 
verwendet werden. Z.B. erlaubt eine quadratische Approximations-Funktion der 
Kostenfunktion die analytische Bestimmung des Minimums der Kostenfunktion in 
jeweiligen Approximations-Bereich. Ahnliche analytische Bestimmungen von Minima 
sind auch mit Cosinus und Sinus Approximationsfunktionen moglich. 

Ausgehend von einem beliebigen Startvektor Xo in dem jeweiligen Approximations- 
Bereich kann das Minimum dann gebildet werden durch die Summe des Startvektors 
x 0 und einer gewichteten Differenz des MeBwerts von der Kalibrierung in Abhangigkeit 
des Startvektors. Konkret kann dies fur das Beispiel eines Luftdatensystems mit 
folgendem Rechenschritt erfolgen: 

x = x 0 + P H R" 1 (x 0 )[y-(c p (x 0 )q co +p so )] 
Zur Erlauterung dieses Losungsschritts wird ausgefuhrt: 

Die Kovarianz Cov(y) einer geeigneten Anzahl von Druckmessungen y als MaBzahl der 
Genauigkeit ergibt sich aus der Losung der Kostenfunktion im Schritt 10: 
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R(x 0 ) = Cov(y) = q c 2 Cov(c p )(x 0 ) + a p 2 

In einem weiteren Schritt ergibt sich die Jacobimatrix der Druckmessungen durch: 

H(x 0 ) = dy/dx(x 0 ) von dem Druckmodell y = q c c p (x) + p s 

Weiterhin ist die Kovarianz bzw. Genauigkeit der Luftdaten x bestimmbar durch: 

P(Xo) = Cov(x) = (H T R- 1 H)" 1 . 

Daraus ergibt sich der angegebene Losungsschritt fur das erfindungsgemaRe 
Verfahren zur Berechnung des lokalen Miniums aus dem Startvektor: 

x = x 0 + P H R" 1 (x 0 )[y-(c p (x 0 )q co +Pso)] 

Urn Rechenzeit zu sparen, kann der Term P H R" 1 (Xq) und die Kalibrierflache c p (x 0 ) in 
Form einer Tabelle abgelegt sein. 
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1 . Verfahren zur Ermittlung einer ZustandsgroBe x aus zumindest einem Sensor-Wert 
mit Hilfe einer fur einen MeBwert y erstellten Kostenfunktion zur Implementierung in 
einer Recheneinheit eines Sensorsystems mit zumindest einem Sensor, wobei die 
Kostenfunktion vom jeweiligen zu messenden Zustand x abhangt und die Abweichung 
eines tatsachlichen MeBwerts von den Kalibrierungen in Abhangigkeit des Zustands x 
angibt, um aus dessen Minimum den gesuchten Zustand x zu ermitteln, 

dadurch gekennzeichnet, 

daB fur die Kostenfunktion zumindest eine Approximations-Funktion auf der Basis 
zumindest eines Approximations-Bereichs innerhalb des Zustandsbereichs x aufgestellt 
wird, durch die eine Approximation der Kostenfunktion mit den Approximations- 
Funktionen mit vernachlassigbarem Fehler erfolgt, wobei die Summe der 
Approximations-Bereiche den gesamten relevanten Zustandsbereich abdecken, 

daB zumindest in einer Auswahl von Appfoximations-Bereichen alle lokalen Minima auf 
der Basis der Approximation bestimmt werden, 

daB gegebenenfalls eine Bestimmung des globalen Minimums aus dem Vergleich der 
lokalen Minima erfolgt. 

2. Verfahren zur Ermittlung einer ZustandsgroBe x aus zumindest einem Sensor-Wert 
nach dem Anspruch 1, dadurch gekennzeichnet, 

daB zur Ermittlung des Minimums der Kostenfunktion eine Bestimmung aller Minima in 
jedem Approximations-Bereich erfolgt. 

3. Verfahren zur Ermittlung einer ZustandsgroBe x aus zumindest einem Sensor- 
Wert nach dem Anspruch 1, dadurch gekennzeichnet, 

daB eine quadratische Approximation der Kostenfunktion erfolgt und zur Ermittlung des 
Minimums der Kostenfunktion in einer Auswahl von entsprechend festgelegten 
Approximations-Bereichen die Minima in jedem Approximations-Bereich zu 
ermoglichen, wobei eine Bestimmung des lokalen Minimums auf der Basis der 
quadratischen Approximationsfunktion erfolgt und jedes lokale Minimum x in der 
Nachbarschaft der ausgewahlten Gitterpunkte ausgehend von einem Startvektor x 0 
bestimmt wird durch die Summe des jeweiligen Startvektors x 0 und einer gewichteten 
Differenz des MeBwerts von der Kalibrierung in Abhangigkeit des Startvektors. 
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4. Verfahren zur Implementierung in einer Recheneinheit eines Sensorsystems nach 
dem Anspruch 1, 2 oder 3, dadurch gekennzeichneit, daB die 
Approximationsfuhktionen und die Nachbarschaften so festgelegt werden, daB die 
Approximationsfunktion in der Nachbarschaft des jeweiligen Gitterpunktes nur ein 
Minimum besitzt, und daB die Bestimmung des lokalen Minimums mittels eines 
Gradientenverfahrens erfolgt. 

5. Verfahren zur Implementierung in einer Recheneinheit eines Sensorsystems nach 
dem Anspruch 1 oder 2, dadurch gekennzeichnet, daB 

die Approximationsfunktion eine Funktion ist, deren Minima mittels analytischer 
Verfahren bestimmbar sind. 



6. Verfahren zur Implementierung in einer Recheneinheit eines Sensorsystems nach 
einem der voranstehenden Anspriiche, dadurch gekennzeichnet, daB 
die zumindest eine Approximations-Funktion und der zumindest eine Approximations- 
Bereich in einem rekursiven Verfahren vor der Aufstellung der Kostenfunktion erfolgt. 



7. Sensorsystem mit zumindest einem Sensor (3) zur Erfassung eines System- 
Zustands und mit einer Recheneinheit (4), wobei die Recheneinheit (4) das Verfahren 
zur Bestimmung einer ZustandsgroBe aus zumindest einem Sensor-Wert nach einem 
der vorabstehenden Anspriiche umfaBt. 
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Zusammenfassung 



Verfahren zur Ermittlung einer ZustandsgroBe x aus zurnindest einem Sensor-Wert mit 
Hilfe einer fur einen MeBwert y erstellten Kostenfunktion zur Implementierung in einer 
Recheneinheit eines Sensorsystems mit zurnindest einem Sensor, wobei die 
Kostenfunktion vom jeweiligen zu messenden Zustand x abhangt und die Abweichung 
eines tatsachlichen MeBwerts von den Kalibrierungen in Abhangigkeit des Zustands x 
angibt, urn aus dessen Minimum den gesuchten Zustand x zu ermitteln, wobei fur die 
Kostenfunktion zurnindest eine Approximations-Funktion auf der Basis zurnindest eines 
Approximations-Bereichs innerhalb des Zustandsbereichs x aufgestellt wird, durch die 
eine Approximation der Kostenfunktion mit den Approximations-Funktionen mit 
vernachlassigbarem Fehler erfolgt, wobei die Summe der Approximations-Bereiche den 
gesamten relevanten Zustandsbereich abdecken, wobei zurnindest in einer Auswahl von 
Approximations-Bereichen alle lokalen Minima auf der Basis der Approximation 
bestimmt werden, wobei gegebenenfalls eine Bestimmung des globalen Minimums aus 
dem Vergleich der lokalen Minima erfolgt. 



(Fig. 1) 



